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Abstract 

By expressing the Holstein-Primakoff transformation in a symmetric form a mod- 
ified spin-polaron technique utihzed on triangular-lattice antiferromagnet is devel- 
oped. With the technique, we have treated an extended t-J model, calculated the 
quasiparticle dispersion, and we also compared the the dispersion with that obtained 
by other method. 
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1 Introduction 



It is found that the cobalt oxide Na2:Co02 ■?/H20 (x ~ 0.35, y ~ 1.3) has a tri- 
angular lattice in the C0O2 planes [1,2, 3, 4]. This material is a fully frustrated 
system when only the nearest-neighbor (NN) correlation is taken into account. 
So it should be necessary to pay more attention on the triangular-lattice anti- 
ferromagnet (TAFM) system. With this motivation, we developed a modified 
spin-polaron technique to discuss the quasiparticle dispersion of the TAFM. 



2 Holstein-Primakoff transformation 



In order to develop a spin-polaron technique on the TAFM, we first express the 
Holstein-Primakoff (HP) transformation in a symmetric form on the square- 
and triangular-lattice, respectively. 
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2.1 Square lattice 



A square- lattice AFM consists two sublattices, one spin-up and the other spin- 
down. We introduce a two-component vector 



(where being boson operators). Then, the HP transformation can be ex- 
pressed in terms of the vector /3j as 



S^^S-alai = S(3}aA = S(3j(^l A, (2) 

=^ (^4 V'2S' - altti + \J2S - alaitt,^ = Spja^pi, (3) 
Sf^'- (^ay2S-alai - yj2S - ala^a^ = -S^a.A, (4) 



or 



si = Sj3lal3,, (5) 
where B is Pauh matrix. 



In spin- wave theory (SWT), in order to introduce only one type Boson, a 
canonical transformation is usually performed to change the Neel configuration 
I till ••• > into a ferromagnetic state with all spins up, i.e., the z axis of spin- 
down sublattice must be upturned, forming the new local coordinate o — x'y'z'. 
Now we investigate how the vector /3j is rotated with the coordinate rotation. 
Suppose the new coordinate is obtained by rotating the old one by 180° about 
its X axis, with z' pointing along the local Neel direction, the direction of 
x'-axis is invariable and y'-axis is pointing along —y. Accordingly, the spin 
components become as 






(6) 
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R 



1 











-1 











-1 



(7) 



2 




is SO (3) matrix. [Sj , Sj , Sj) are spin components in the old coordinate, and 
{S'j', S'j , S'f) in the new local coordinate. 

With the coordinate rotation, /dj become Pj. We suppose that 1) they are 
related through a indeterminate matrix u{R): 

= u{R)f3j, (8) 
and 2) the HP transformation is unchanged in its form, i. e., 

^^^S^^a^^. (9) 
Then, we have immediately the relation 

I S(3]u\R)a,u{R)(3,\ 
= Sl3]n){R)ayu{R)l3j . (10) 
\S(5]u\R)a,u{R)(5j 

Prom this equation the indeterminate matrix u{R) can be easily solved, and 
the results is 



<R) = {\ I)- (11) 

Because the new coordinate is fixed on the spin-down sublattice and the old 
one on the spin-up sublattice, the vector /3j has the form of Eq. (1) on spin-up 
sublattice, and the form 



/*;^«m^-i=(^^^) (12) 

on spin-down sublattice. If the prime is omitted and the spin-up and -down 
sublattices are distinguished by indices, we have 



_^ ( ^25 g gpi^ _ sublattice), 

A= V / (13) 

( / ' — I — I (i e spin — down sublattice). 

The HP transformation can be merged into an unison form on the both sub- 
lattices: 
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(14) 



with a = i,j corresponding to spin-up and -down sublattices, respectively. It 
is easily verified that on both sublattices the two-component vector satisfies 
the normal condition 

PlPi = 1- (15) 

2.2 Triangular lattice 



Analog on the square-lattice AFM, now we express the HP transformation 
on the TAFM. Unlike the square-lattice AFM, the TAFM has three sublat- 
tices(called A, B and C) with three 120°-Neel states, and their local coordi- 
nates can't be simply divided into spin-up and spin-down sublattices, but into 
three. 

Following Miyake [5,6], we define the local (spatially varying) coordinates o- 
x'y'z', with y' pointing along the old z direction and z' pointing along the 
local 120°-Neel direction. When x' is rotated by 0°, 120° and 240° about y' (z) 
axis respectively, three new local coordinates are formed, which are fixed on 
the sublattices A, B and C, respectively. In the three new coordinates a spin 
operator has three forms: 

r {Sf, St, (i e A) 

{ST. S'n = i-fst - |5f , St, -^St + fsf) {i e B) (16) 

(Vlqx_iqy qz _iqx _ Viqy\ r<\ 
We merge the three form into one 



lsr\ 




( s^\ 




= K' 


sA {a = A,B,C), 









(17) 



Then the matrix R^^ can be easily resolved from the Eqs.(16), and the inverse 
matrices are 










■ 







Vo 


1 





Ra=\1 Q Q\,Rb=\ -\ ^ 



2 2\ 12^ 2' 

(18) 



2 ^ 2 » 



/ 



,Rc= I -1 

3 1 / 



Similarly in the Section 2.2, we here also introduce a two-component vector 
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in the old coordinate, and the HP transformation is still expressed in terms 
of Pi{0) as 



Si^S(3liO)a(3iiO). (20) 

— * 

When the coordinate is rotated, spin operator changes from si to s'^, and the 
introduced matrix from (3i{0) to Pi{a) (where a — A, B C corresponding to 
the three new coordinates, respectively). Wc suppose the HP transformations 
on the new coordinates are expressed in an unison form 

^,^S(3}(a)a(3,(a), (21) 
And we suppose also that 

Pi{a) = M(i?„)A(0). (22) 
Prom the Eqs.(21) and (22), we can express s'i in terms of /3i(0): 



s'ia = SPj{0)u\R^)au{Ra)Pi{0) {a G A, B, C). (23) 

The indeterminate matrix u{Ra) can be determined by substituting the Eqs. 
(20) and (23) into (17), and the results is 

. „ \ { cos a — sin q; \ 
uiRa) = . (24) 
^ ^ Vsmo; cosa / ^ ' 

with a = 0,27r/3, — 27r/3 on sublattices A, B and C, respectively. Eventually, 
Eq. (21) is just the HP transformation in the three local coordinates of TAFM. 



3 Modified spin-poleiron technique 



After expressing the HP transformation in terms of the introduced matrix, we 
now develop a modified spin-polaron technique. The spin-polaron picture was 
proposed early by Schmitt-Rink Varma and Ruckenstein [7] to deal with the t- 
J model on square latticc[8,9,10,ll]. In this picture the electron-annihilation 
operators are expressed as pure hole operators or composite operators, for 
example. 
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(25) 



with s\ being the hard-core Bose operators. A similar spin-polaron picture was 
proposed by Liu and Manousakis [8] by introducing two types of holes and 
two types of spinous on spin-up and spin-down sublattices, respectively. 

Since the electronic operators Qo- {Cl„) appear always in pairs in physical 
quantities (for example, the kinetic operator '}Z<ij>,aCia^j<yi current opera- 
tor '}2<ij>,a RiClcrCjai the Hamiltoniau if, and for the t-J model, the single- 
occupancy constraint Z^o- C'lo-Qo- < !)> we should deal directly with the pair 
operator Y,^ Cj^Cja — Cj-^Cj-i + Cj^Cji, rather than the single electronic oper- 
ators Ci^ (ct). 

Because the electron hopping operators Cj-^Cji and Cj^Cji correspond to the 
same hole hoppings from the site i to j, the term Cj-^Cj^ + Cj^Cji should be 
proportional to the hole hopping operators hihp or 

^iT^jT + ^ii^ji = i^ijhih]. (26) 

The factor Kij should be related to boson operators ai(j) and aj^^-^, and one 
may expand it in terms of a series of these boson operators, 

Kij = Ao Ai{a\ ttj) A2{a\aj ajcj) (27) 

where Aq,Ai,... are indeterminate coefficients. Determination of them is de- 
termination of the modified spin-polaron technique. 

On the one hand, in terms of the electron operators and the Pauli matri- 
ces, the spin operators can be expressed as = \ '}Zaa' ^la^aa'Cia' , and the 
corresponding 2;-component reads 

^ ^ia^^aa'Cia' = S{C\^Ci^ - C/^C^). (28) 

aa' 

On the other hand, the component can be expressed as 
So we have the relation 
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If the negative sign is changed as positive, one immediately has 



(4QT + 4Ca)=/3l(J (29) 

This is exactly true as it is an identity. Enlightened by this relation, we may 
extend it from the same site to different site: 



iChC,^ + ClC,^)^f^}[l ^^)Pj = PlPj- (30) 
This extension implies that the factor Kij have been selected as 



and the coefficients as Aq — 0, Ai = -^,^2 = 0, ^3 = • • •. Finally, the Eqs. 



(26) and (31) make up the modified spin-polaron transformation. 

It should be stressed that there may be other selections to A's. For example, 
one may suppose Kij — f{Pi,Pi,Pj,(3j) as long as the operators Kij satisfy 
the necessary requirements such as conjugation for the permutation of i and 
j, and unitarity when i = j. Different selection may correspond to different 
magnon-holon coupling strength. 

Now we rewrite the modified spin-polaron transformation in a compact form 



Y,ClC,, = h4{a)h]Pj{P), (32) 

(7 

where the index a {P) is for distinguishing different sublattices with the site i 
(j) belonging to the sublattice a (P). 

Because A (a) satisfies the normal condition 



f3j{a)f3i{a) = 1, (33) 

on the same site, with the modified spin-polaron technique the no-double 
occupancy constraint is automatically built in: 



^ Cld, = hiPl{a)hlPi{a) = hih\ < 1. (34) 
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4 Application of the spin-polaron technique on TAFM 



Now we use the modified spin-polaron technique to treat the TAFM. Here we 
use the extended t-J model to describe its physics. Then, when the long-range 
correlations are taken into account, the Hamiltonian reads 



H — Htt' + Hj, 

Htt'——t ^ C^^Cjcr — t' ^ C^^Cja — fi^C^^Cia, (35) 

{ij)icr {ij)2T i 

Hj^ J J2S,-Sj, (36) 

where the summations and {i,j)2 run over the NN and next-nearest 

neighbor (NNN) pairs respectively and the operators Cj^ are subjected to the 
single-occupancy constraint. 

The spin-spin correlation part Hj of the Hamiltonian can be treated with the 
HP transformation. In k space the free part of the spinon energy is 



Hj = J2'^kaiak, (37) 

k 

where ccfe are spinon operators. The spin- wave dispersion is 



= ^W[(l + 27i'^)(l-7i'^)] (38) 

where 7^^^^ = | T^sW e'^'^^^^ is the summations over the NN sites. And the vec- 
tors covers the six NN neighbors etr, — etr, —^Cx + ^Sy, \ex—'^ey, —\ex — 

^e.y and |e-r -|- ^e^,), e.^ being one of the basis vectors, and normal to e^- 
Eq. (38) is exactly the same as that obtained by Leung and Runge [6]. 

With the transformation Eq. (32) the Hamiltonian can be expressed by 
boson and hopping operators. If we preserve the second order of bosons, the 
it reads 



Htt' —Ht + Hf, 

J2 hih]{al-aj)- ^ hihliaj - aj) 

— -^t X! hihUa-ai + a]aj — 2ajaj) — A* hihl + H.c, (39) 
8<b / • v 
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H't ~ -t' J2 



t 
j 

{y>2 



+ H.c. (40) 



In k space with Bogliubov transformation, we have 

Htt = J2^khihk + H', (41) 

k 

where hk are holon operators. The first term describes the holon hopping, and 
holon dispersion is 



^fc--^[^7i'^-2t'7f], (42) 

where 7^^'* = - J2s(2) e'^'^'^''^^ is the summations over the(NNN) sites. In Eq. (41) 
the second term H' describes the interaction between the holons and spinons. 



H'^E {Vi,h,hlal_^ + V,,h,hlak^,) , (43) 

kp 

where Vkp is the coherence factors. Here we will not discuss it in detail, but 
pay attention mainly on the holon dispersion. 

Eq. (42) gives out the holon dispersion when both NN and NNN hoppings are 
included. If the NNN hopping is ignored, the spectrum reduces to 



It is a periodical function. Its amplitude is one half of Trumper's [12] and only 
one sixth of Azzouz's [13]. This means that the present dispersion is the least. 
Why? We know that the TAFM is fully frustrated, and the ground state is 
very disordered. The disorder certainly flattens the dispersion. So the property 
of spin frustration is more fully maintained within the present theory. 

In summary, after introducing a two-components matrix we express the HP 
transformation in a symmetric form, based on this we developed a modified 
spin-polaron technique. With the technique we calculated the quasiparticle 
dispersion of an extended t-J model. The dispersion is more reasonable than 
that obtained by other methods. The present theory can fully describe the 
frustrated TAFM. 
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